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Abstract — In this paper we develop a new statistical model for 
the irradiance fluctuations of an unbounded optical wavefront 
(plane and spherical waves) propagating through a turbulent 
medium under all irradiance fluctuation conditions in homoge- 
neous, isotropic turbulence. The major advantage of the model is 
that leads to closed-form and mathematically-tractable expres- 
sions for the fundamental channel statistics of an unbounded 
optical wavefront under all turbulent regimes. Furthermore, it 
unifies most of the proposed statistical models for the irradiance 
fluctuations derived in the bibliography providing, in addition, 
an excellent agreement with the experimental data. 

Index Terms — Atmospheric turbulence, free-space optical com- 
munications, optical intensity modulation, scintillation. 



I. Introduction 

Atmospheric optical communication is receiving consider- 
able attention recently for use in high data rate wireless links. 
However, even in clear sky conditions, wireless optical links 
experience fluctuations in both the intensity and the phase of 
an optical wave propagating through this medium H i due to 
time varying inhomogeneities in the refractive index of the 
atmosphere. The reliability of an optical system operating in 
such environment can be deduced from a mathematical model 
for the probability density function (pdf) of the randomly fad- 
ing irradiance signal (scintillation). For that reason, one of the 
goals in studying optical wave propagation through turbulence 
is the identification of a tractable pdf of the irradiance under 
all intensity fluctuation regimes. 

Over the years, many irradiance pdf models have been 
proposed with different degrees of success. Perhaps the most 
successful models are the log-normal [2], the lognormal- 
Rician (or Beckmann model) (|3] and the gamma-gamma H. 
The scope of the lognormal model is restricted under weak 
irradiance fluctuations. The Beckmann model, although shows 
excellent agreement with experimental data, however a closed- 
form solution for its integral is unknown in addition to its 
inherent poor convergence properties. Finally, the gamma- 
gamma pdf was suggested as a reasonable alternative to Beck- 
mann's pdf because of its much more tractable mathematical 
model. 

Now, through this paper, we propose a new and generic 
propagation model and, from it, we obtain a new and unifying 
statistical model for the irradiance fluctuations. The proposed 
model is valid under all range of turbulence conditions (weak 
to strong) and it is found to provide an excellent fit to 
experimental data. Furthermore, the statistical model presented 
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here can be written in a closed-form expression and it contains 
most of the statistical models for the irradiance fluctuations 
that have been proposed in the bibliography. Thus, the utility of 
our model can be extended to many other areas of application 
as the study of fading in radiofrequency channels. 



II. The proposed model of propagation 

Assume an electromagnetic wave is propagating through 
a turbulent atmosphere with a random refractive index. As 
the wave passes through this medium, part of the energy is 
scattered and the form of the irradiance probability distribution 
is determined by the type of scattering involved. In the physical 
model we present in this paper, the observed field at the 
receiver is supposed to consist of three terms: the first one 
is the line-of-sight (LOS) contribution, Ul, the second one 
is the component which is quasi-forward scattered by the 
eddies on the propagation axis, Ug, and coupled to the LOS 
contribution; whereas the third term, U^, is due to energy 
which is scattered to the receiver by off-axis eddies, this 
latter contribution being statistically independent from the 
previous two other terms. The inclusion of this coupled to 
the LOS scattering component is the main novelty of the 
model and it can be justified by the high directivity and 
the narrow beamwidths of laser beams in atmospheric optical 
communications. The model description is depicted in Fig. [T] 

Accordingly, we can write the total observed field as: 
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U ={Ul + U^ expix + jS) 



where Ul = VcV^Gxp (jtfiA), 
and t/G=^(l 1 p)U'g, being 
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independent stationary random processes. Of course, Ul and 
Ug are also independent random processes. In ([T), G is a 
real variable following a gamma distribution with £'[G]=1. 
It represents the slow fluctuation of the LOS component. 
Following |5|, the parameter i7=£'[|[/ip] represents the 
average power of the LOS term whereas the average 
power of the total scatter components is denoted by 2bo = 
E[\Ug\'^+\Ug\'^]. (I) A and 4>b are the deterministic phases of 
the LOS and the coupled-to-LOS scatter terms, respectively. 
On another note, 0<p<l is the factor expressing the amount of 
scattering power coupled to the LOS component. Finally, Ug 
is a circular Gaussian complex random variable, and x and S 
are real random variables representing the log-amplitude and 
phase perturbation of the field induced by the atmospheric 
turbulence, respectively. A plausible justification for the 
coupled-to-LOS scattering component, Ug , is given in 0. 
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Fig. 1. Proposed propagation geometry for a laser beam. 

III. Malaga (A^) probability density function 

From Eq. the observed irradiance of the proposed 
propagation model can be expressed as: 

2 
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Ul 



exp (2x) = YX, 

2 

(small-scale fluctuations) 



X = cxp (2x) 



(large-scale fluctuations), 



where the small-scale fluctuations denotes the contributions to 
scintillation associated with turbulent cells smaller than either 
the first Fresnel zone or the transverse spatial coherence radius, 
whichever is smallest. In contrast, large-scale fluctuations of 
the irradiance are generated by turbulent cells larger than that 
of either the Fresnel zone or the so-called "scattering disk", 
whichever is largest. From Eq. ([l}, we can rewrite the lowpass- 
equivalent complex envelope as: 



R{t) =(c/l 



U§ 



U§ 



Qexp{j4}A) 



+ VPV^exp(j?!)B) + v/(l - p)U's 



(3) 



SO that we have the identical shadowed Rice single model em- 
ployed in |5 1, composed by the sum of a Rayleigh random pha- 
sor (the independent scatter component, Ug) and a Nakagami 
distribution, \/G. Then, we can apply the same procedure 
exposed in |5 | consisting in calculating the expectation of the 
Rayleigh component with respect to the Nakagami distribution 
and then deriving the pdf of the instantaneous power Hence, 
the pdf of Y is given by: 
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(4) 



where /3 — (i?[G])^/Var[G] is the amount of fading parameter 
with Var[-] as the variance operator, being | • | the absolute 
value. We have denoted 7 — 2bo{l — p) and 51 = il + p2bQ + 
2\/2boflpcos {(pA — 4>b)- Finally, iFi (a; c; x) is the Kummer 
confluent hypergeometric function of the first kind. 

Otherwise, the large-scale fluctuations, X = cxp (2%), is 
widely accepted to be a lognormal amplitude [3| but, as in 
H, Q, this distribution is approximated by a gamma one, 
this latter with a more favorable analytical structure. Then: 



fx{x) 



r(a) 



exp {—ax) , 



(5) 



where a is a positive parameter related to the effective number 
of large-scale cells of the scattering process, as in 2]. 



Now, the statistical characterization of / can be obtained. 
It is said that / follows a M distribution if X and Y are 
random variable distributions according to Eqs. (|5) and (2), 
respectively, with /3 being a natural number (a generalized 
expression for /? being a real number can be also derived, 
with an infinite summation, but it is less interesting due to the 
high degree of freedom of the proposed distribution). Then, 
the pdf is represented by: 



k=l \ 
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In Eq. (|6]l, Ky{-) is the modified Bessel function of the 
second kind and order v. In the interest of clarity, the algebraic 
manipulation to prove this result is moved to Appendix A. 

Finally, the fc*'' centered moments of the irradiance signal 
following a M. distribution, denoted by ruk {I), are given by: 
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The proof of this result is, again, moved to Appendix B. 

To conclude this section, we attach in Table H] a list of most 
of the existing distribution models for atmospheric optical 
communications and how they can be generated from the 
distribution model proposed here. The analytical expression 
for the different distribution models include in Table |I] can be 
consulted in the references provided next to their names. 

TABLE I 

Generation of different statistical models froivi the proposed 
m distribution model. 



Distribution model 


Generation 


Rice-Nakagami | r| 


p = 




Var[|C/L|] =0 


Lognormal (T], (7] 


p = 


Var[|C/L|] =0 




7^0 


Gamma 14] 


p = 




7 = 


Shadowed-Rician distiibution 15J 


Var[|X|] = 


K distribution|41, (Sj 


n = and p = 




ox 13 = 1 




n = 


HK (Homodyned K) distribution 8| 


p = 




X = 7 




n = 


Exponential distribution (4] 


p = 




a — >■ 00 


Gamma-gamma | 4 | 


p = 1, then 7 = 


distribution 


q! = 1 


Gamma-Rician distiibution |3 


/3 — >■ 00 
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Fig. 2. The pdf of the scaled log-irradiance for a plane wave (Figures 
(a), (b) and (c)) and a spherical wave (Figure (d)) in the case of: (a) weak 
irradiance fluctuations (af =0.1 and lo/Rp = 0.5); (b) moderate irradiance 
fluctuations (af = 2 and lo/Rp = 0.5); (c) strong irradiance fluctuations 
(it? = 25 and lo/Rp = 0.5); and (d) strong irradiance fluctuations 
^'^Rytov ~ ^ ^^'^ 1(1 /Rf = 0.5). The blue open circles represent simulation 
data, the dashed red hne is from the gamma-gamma pdf with a and /3 
predicted in |9| and the solid black line is from our M distribution model. 

IV. Comparison with experimental data 

In this section, we compare the distribution model 
with some of the published numerical simulation data plots 
in (|9] and ifTOl of the log-irradiance pdf, covering a range 
of conditions that extends from weak irradiance fluctuations 
far into the saturation regime. It is remarkable the evident 
analytical tractability by directly employing Eq. (|6}, with a 
finite summation of /3 terms. Some results are shown in Fig. 
ID emphasizing its extremely high accuracy. Thus, we plot in 
Figs. |2](a)-(c) the predicted log-irradiance pdf associated with 
the M distribution (black solid line) for comparison with 
some of the simulation data illustrated in Figs. 4, 5 and 7 
of ||9]. The simulation pdf values are plotted as a function of 
(In J— < In/ >)/(T, as in f^l, where < In/ > is the mean 
value of the log-irradiance and a — V'^hi/' '■^^ latter being the 
root mean square (rms) value of In /. For sake of brevity, and 
as representative of typical atmospheric propagation, we only 
use the inner scale value = O.bRp so we can include the 
effect of Iq in our results, where the quantity Rp is the scale 
size of the Fresnel zone. We also plot the gamma-gamma pdf 
(red dashed line) obtained in for the sake of comparison. In 
Fig.|2](a) we use a Rytov variance — 0.1 corresponding to 
weak irradiance fluctuations, in Fig. |2](b) we employ — 2 
corresponding to a regime of moderate irradiance fluctuations 
whereas in Fig.|2](c), was established to 25 for a particular 
case of strong irradiance fluctuations. Values of the scaling 
parameter a required in the plots for the M pdf are obtained 
from Andrews' development flU in the presence of inner scale. 

Finally, in Fig.|2](d) we have obtained a very good fitting to 
the simulation data for a spherical wave in the case of strong 



irradiance fluctuations. Following Hill's representation ifTOl . 
the simulation pdf data and pdf values predicted by the Ai 
distribution are displayed as a function of (In / + 0.5(T^)/(t, 
where a was defined above. In this particular case of prop- 
agating a spherical wave, an additional parameter is needed: 
the Rytov parameter, crjzjytov^ defined as the weak fluctuation 
scintillation index in the presence of a finite inner scale. 

For the particular case displayed in Fig. |2](d), the gamma- 
gamma pdf does not fit with the simulation data and, even 
more, the Beckmann pdf did not lend itself directly to numer- 
ical calculations and so are omitted. Nevertheless, the M pdf 
shows very good agreement with the data once again, with 
the advantage of a simple functional form, emphasized by the 
fact that its (3 parameter is a natural number, which leads to 
a closed-form representation. 

V. Concluding remarks 

In this paper, a novel statistical model for atmospheric 
optical scintillation is presented. Unlike other models, our 
proposal appears to be applicable for plane and spherical 
waves under all conditions of turbulence from weak to super 
strong in the saturation regime. The proposed model unifies 
in a closed-form expression the existing models suggested 
in the bibliography for atmospheric optical communications 
(see Table HJi. In addition to the mathematical expressions and 
developments, we have introduced a different perturb ational 
propagation model that gives a physical sense to such existing 
models. Hence, the received optical intensity is due to three 
different contributions: first, a line of sight component, second, 
a coupled-to-LOS scattering component, as a remarkable nov- 
elty in the model, that includes the fraction of power traveling 
very closed to the line of sight, and eventually suffering from 
almost the same random refractive index variations than the 
LOS component; and third, the scattering component affected 
by refractive index fluctuations completely different to the 
other two components. The first two components are governed 
by a gamma distribution whereas the scattering component is 
depending on a circular Gaussian complex random variable. 
All of them let us model the amplitude of the irradiance 
(small-scale fluctuations), while the multiplicative perturbation 
that represents the large-scale fluctuations is approximated to 
follow a gamma distribution. 

Finally, we have made a number of comparisons with 
published plane wave and spherical wave simulation data over 
a wide range of turbulence conditions (weak to strong) that 
includes inner scale effects. The distribution model is found 
to provide an excellent fit to the simulation data in all cases 
tested making computations extremely easy due to its closed- 
form expression. 

Appendix A: proof of Equation dU 

To prove Eq. (|6]l, we can obtain the Laplace transform, 
£[/y(?/);s], of the shadowed Rice single pdf, fviy), writ- 
ten in Eq. (IHi, in a direct way, with the help of Eq. (7) 
of Ref. 15], since the moment generating function (MGF) 
and the Laplace transform of the pdf friy) are related by 
M[fY{y);-s]^C[fY{y);s]: 
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(9) 



Now, let us consider the following Laplace-transform pair 

r (ly+l) (s - A)" (s - fiy^^ ^ n!t''-"e'^*Lr" [(A - m) ^] , 

Re(i^) > n — 1; 

(10) 

given in ifTTI . Eq. (4) in pp. 238, where the minor error in the 
sign of the argument of the Laguerre polynomial found and 
corrected in [12] has already taken into account. If we denote 

n = /3 — 1 and 1/^/3—1, then 

/3-1 / a \ -P 



^ 7' 7/3+Jl' 



(/3-1)! (s + i 



/3 



7/3 + 



n t 



(11) 



where is the Laguerre polynomial of order n. If we 

substitute Eq. (fTTl i into Eq. then the pdf of Y can be 
expressed as: 



fy (y) = - 



7/3 



7 V7/3 + n' 



Gxp 



/3 



7/3 + 
(7/3 + 7 



(12) 



Now, to obtain the unconditional Ai distribution, from Eqs. 
© and (fT2l l. we can form: 

//(/) - / fY{I\x)fx{x)dx 
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exp 
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7/3 + f7' a: 
exp (—ax) dx. 



(13) 



7/3 + SI' 7x 

By expressing the Laguerre polynomial in a series. 



(14) 



as was shown in Eq. (8.970-1) of Ref. 113.1 . and using Eq. 
(3.471-9) of Ref. Qll, it follows that Eq. O becomes 



fiil) 



r(a)7i+t \-fl3 + n' 
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Z-^^\k - 1 J k - v. l7(7/3 + r2') 
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Finally, Eq. ( fTsT i can be rewritten as Eq. (|6), as was already 
anticipated in Section III. 

Appendix B: proof of Equation ^ 

From Eq. (jJJ, the observed irradiance, /, of our proposed 
propagation model can be expressed as: / = XY, where the 
pdf of variables X and Y were written in Eqs. ^ and (|4]l, 
respectively. Based on assumptions of statistical independence 
for the underlying random processes, X and Y, then: 

mfc (/) = E [X'^] E = ruk {X) (F) . (16) 

From Eq. (2.23) of Ref. lfT4l . the moment of a Nakagami-m 
pdf is given by: 

r(a + fc) 



rrik (X) 



(17) 



r (a) OL^ 

On the other hand, from Eqs. (fTST i and (fl4] i. and employing 
Eq. (3.381-4) of Ref. QSl, obtain the moment of the Rician- 
shadowed distribution, given by: 



ruk [Y) = 



7/3 



1 



7 V7/3 + f^' 
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r=0 



/3-n 1 



r(fc-f-r + l) 



(18) 
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Finally, when performing the product of Eq. ([TtI by Eq. ([18), 
we certainly obtain Eq. dHJ. 
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